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ABSTRACT 

The nuclear reactions that occur in the stellar progenitors of white dwarfs (WDs) lead to an internal composition 
of 12 C, 16 0, and a "contaminant" nucleus, 22 Ne . The 22 Ne is produced by helium captures on 14 N left from 
hydrogen burning via the CNO cycle. By virtue of its two excess neutrons (relative to the predominant A = 2Z 
nuclei), a downward force of rj 2m p g is exerted on 22 Ne in the WD interior. This biases its diffusive equilibrium, 
forcing 22 Ne to settle towards the center of the WD. We discuss the physics of the gravitational settling when the 
WD is in the liquid state and the luminosity generated by it. This modifies the cooling of WD's with masses in 
excess of M . The current uncertainties in the microphysics even allow for solutions where a 1 .2M WD remains 
mostly liquid for a few Gyrs due to the internal heating from 22 Ne sedimentation. This highlights the need for 
an accurate calculation of the inter-diffusion coefficient, especially in the quantum liquid regime relevant for high 
mass WD's. There is also time in old, liquid WD's (such as those found in cataclysmic variables and possibly in 
accreting Type la progenitors) for partial settling. 

Subject headings: diffusion - novae, cataclysmic variables - stars: abundances, interiors - supernovae: general - 
white dwarfs 



1. INTRODUCTION 

The slowest step in the CNO cycle is the proton capture on 
14 N. This results in all CNO catalysts piling up at 14 N when 
H burning is completed during the main sequence. During He- 
lium burning, the reactions I4 N(a,7) 18 F(/? + ) 18 0(a,7) 22 Ne en- 
sue and convert the 14 N into 22 Ne , resulting in a mass frac- 
tion of 22 Ne , X22 ~ Zcno ~ 0.02 for Population I. The pos- 
sibility that this nucleus plays an important role in the ener- 
getics of cooling WD's was first noted by Isern et al. (1991). 
They discussed the possibility of gravitational energy release 
if 22 Ne phase separates at crystallization. This was followed 
by more detailed studies (Xu & Van Horn 1992; Ogata et. al. 
1993; Segretain et al. 1994; Segretain 1996) which differed in 
their conclusions about the final state once crystallization was 
complete; ranging from all the 22 Ne in the stellar center to most 
of it having an unchanged profile. 

What about 22 Ne sedimentation in the liquid state? Ogata 
et al. (1993) showed that demixing in the liquid state never 
occurs for the C/0/ 22 Ne liquid (namely it will remain in the 
fully mixed state when the temperature is safely above melting). 
This allows us to focus on the gravitational settling of individ- 
ual 22 Ne nuclei in the liquid parts of the white dwarf (Bravo et. 
al. 1992). This happens prior to the onset of cystallization in 
a cooling WD and during the lifetime of a WD that remains in 
the liquid state due to heating from accretion. 

Degenerate electrons dominate the pressure in the WD in- 
terior, so that the upward pointing electric field at mass co- 
ordinate m{r) is eE w 2m p g, where g = Gm(r) jr 2 and m p is 
the proton mass. The downward force on a 22 Ne nucleus 
is F = llnipg- \0eE = 2m p g, giving a potential energy drop 
U(r) = J^2m p gdr between radius r and the center. Across 



the whole star, this amounts to 285 keV for M = O.6M and 
1600 keV for M = 1.2M . It is the excess neutrons of the 
22 Ne nucleus (relative to A = 2Z) that make it special (the same 
is true of the less abundant 56 Fe). Differentiating forces on the 
carbon and oxygen nuclei only arise from the slight differences 
in masses from Am p due to the nuclear binding, which is orders 
of magnitude smaller than the force on 22 Ne . 

The potential energy (the energy that would be released if 
all the 22 Ne sinks to the center) of a uniformly mixed model 
(X 22 = 0.02 is constant) is 1.4 x 10 47 ,3.2 x 10 47 ,6.8 x 10 47 and 
1 .5 x 10 48 erg for masses of 0.6, 0.8, 1 .0, &1.2M . For 0.6 and 
1.0M o , this corresponds to 5.8 x 10 15 & 1.7 x 10 16 erg gr" 1 of 
possible energy release from 22 Ne settling. At a temperature of 
10 7 K, the specific heat per oxygen ion is 7.7 x 10 13 erg gr" 1 , 
so that the energy stored in 22 Ne is comparable to the thermal 
energy content of a cooling WD at late times. 

In §2, we estimate the inter-diffusion coefficient for a trace 
22 Ne nucleus in the liquid C/O interior. The settling timescale 
for 22 Ne is ~ 10 Gyr, making it critical for a more accurate 
calculation of the inter-diffusion coefficient. The settling time 
decreases for more massive WDs. In §3, we speculate on the 
modifications the 22 Ne settling is likely to have on accreting 
WD's and maybe Type la supernovae. We calculate the lumi- 
nosity from the gravitational energy release of sinking 22 Ne in 
§4. This affects the cooling of WD's prior to crystallization and 
can possibly heat them to a level hot enough to maintain them 
in the liquid state for a Hubble time. 

2. DIFFUSION AND SETTLING OF THE 22 Ne 

We presume that the WD interior consists of a single element 
(either 12 C or 16 0) of mass Am p and charge Ze, and that 22 Ne is 
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a trace component at density «n 6 = X22p/22m p . Our WD is con- 
structed purely from the degenerate electron equation of state. 
For the properties we need for the 22 Ne diffusion problem, it 
is a fine approximation to neglect the ionic contribution to the 
pressure. 

Relative diffusion of charged species in a correlated liquid 
is substantially different than the diffusion problem in the WD 
atmosphere. In those less dense regions, traditional Coulomb 
cross sections and ideal gas equations of state dominate (e.g. 
Fontaine & Michaud 1979; Vauclair & Vauclair 1982). As dis- 
cussed by Paquette et al. (1986), there is substantial uncertainty 
in the diffusion coefficients in the regions of the WD where 
the ion Coulomb coupling becomes strong. For a classical one 
component plasma (OCP) with ion separation, a, defined by 
a 3 = 3/47171;, where n,- = p/Am p , the importance of Coulomb 
physics for the ions is measured by 



(Ze) 2 1/3 



10 7 K 



1/3 
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Since kT <C U(r) over most of the stellar radius (equivalent 
to saying that the 22 Ne scale height is much less than the WD 
radius), most of the drift of the 22 Ne is dominated by "falling" at 
the local speed set by the diffusion coefficient V ks 2m p gD/kT, 
which with equation (||) gives 



V = 



18m p g 
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In mass coordinates, m(r), the time it takes to fall from m\ to 

m.2 is 

Zer 1 / 3 dm(r) 
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For a constant density star, this becomes At w t s J dlnm(r), 
where 
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where /% = p/10 6 grcm" 3 . Crystallization occurs when V ex- 
ceeds 173 (see Farouki & Hamaguchi 1993 and references 
therein). We only consider settling of 22 Ne in the liquid state, or 
1 < r < 173, the presumption being that the dramatic increase 
in viscosity expected in the solid state will prohibit further grav- 
itational sedimentation. 

The inter-diffusion coefficient, D, of 22 Ne in a C/O mixture 
with r ~ 1 - 100 has not been specifically calculated and so we 
construct here our best estimate. One way to estimate D is to 
use the Stokes-Einstein relation for a particle of radius a p (in 
this case 22 Ne ) undergoing Brownian motion in a fluid of vis- 
cosity f], which gives D = kT j\-na p r] when the fluid is allowed 
to "slip" at the particle/fluid interface. This estimate works well 
in other liquids at small scales, where it is accurate at atomic di- 
mensions (Hansen & McDonald 1986). An alternate route is to 
use the self-diffusion coefficient in the OCP 



D « 3w„fl 2 r- 4/3 , 



(2) 



Using the central densities as an estimate, a pure oxygen 0.6M Q 
WD at T = 10 8 K (this temperature is appropriate for rapidly ac- 
creting WD's) has t s sa 19.6 Gyr, so that complete sedimenta- 
tion is not likely for the most common 0.6M Q WD's. How- 
ever, the strong density dependence means that this settling 
time decreases for massive WD's. For example, a pure oxygen 
1.2(1.3)M WD at 10 8 K has f s « 4.8(3)Gyr, in rough agree- 
ment with the previous simple estimates of Bravo et al. (1992). 

Figure 1 shows the time it takes for a 22 Ne nucleus that starts 
at the surface (m(r) = M) to fall to the location, m(r). This 
is an integration of equation (Q) with a model WD. The solid 
(dashed) lines are for pure oxygen (carbon) and bracket the C/O 
WD's. We show two masses, 0.6 and 1.2 M G , making it clear 
that 22 Ne settles faster in a massive WD. 

3. ACCRETING WHITE DWARFS AND TYPE IA PROGENITORS 

Unlike cooling WDs, the interiors of massive, accreting 
WD's are maintained in the liquid state by the heating from 
accretion. For slowly accreting (M ~ 10~ 10 Moyr -1 ) WD's in 
cataclysmic variables below the period gap, the core tempera- 
tures are in the range T « (1— 3) X 10 7 K (Nomoto 1982), so 
that they remain liquid while accreting for a few Gyrs. A 1M Q 
pure oxygen WD in such a setting will undergo neon settling in 
the outer layers as t s « 12.6 Gyr at T = 3 x 10 7 K. 

Though the debate rages regarding appropriate Type la pro- 
genitor models, most favor rapidly accreting WDs that ignite 
carbon once the mass is near the Chandrasekhar value. The ac- 
cretion rates considered are too rapid to allow 22 Ne settling dur- 
ing this phase of accretion. However, prior evolution at lower 
accretion rates or an extended period as a cooling WD will have 
allowed for some Neon settling. Though unclear just how im- 
portant this will be, the effects of 22 Ne settling for the Type 
la ignition of a near-Chandrasekhar mass WD are important to 
note. The most obvious one is that 22 Ne settling will develop 
a gradient in the electron mean molecular weight. This would 
modify the convective criterion in the ignited core and possibly 
affect the evolution from the thermal runaway to a propagating 
flame front (Garcia-Senz & Woosley 1995). 1 

'Our present understanding of the microphysics does not predict the development of a pure 22 Ne core. However, if this were to occur, it is unlikely to modify 
the ignition conditions, as the electron energy threshold for the reaction 22 Ne +e~ — > 22 F+u e is 1 1.36 MeV, which requires a density in excess of 2.1 X 10 10 g cm . 
Carbon ignites at densities about an order of magnitude less (Nomoto, Thielemann & Yokoi 1984), so that the initial trigger of carbon burning is unchanged. However, 
a pure 22 Ne core with a mass of 0.028Mq would have a radius of R^ e f» 190 km at a density of 2 X 10' g cnr 3 , forcing the C/O mixture to ignite in the shell covering 
the 22 Ne core. 



calculated by Hansen, McDonald and Pollock (1975), where 
2 _ //(,„„ ;o tu "i™™- fm„,m„„, How 
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uj p = AimiiZe) /Am p is the ion plasma frequency, 
well does this agree with the Stokes-Einstein estimate? 
ing the numerically derived viscosity of the OCP, r\ rj 
Q.lpu p a 2 (T/lQ) 03 (Donko & Nyiri 2000 and references 
therein) and setting a = a p , the derived D = kT /Airarj agrees 
with equation (||) to about 20% for all V > 10. 

We also compared equation to Fontaine's (1987) rough 
fitting formula for the interdiffusion coefficients and found 
agreement to w 30%. So, for now, we will use equation (^|) for 
our estimate of D, keeping in mind that there has yet to be an 
appropriate microphysical calculation. We will show later that 
there is astrophysical motivation for a more accurate calculation 
of the inter-diffusion coefficients, especially in the quantum liq- 
uid regime that is relevant for the massive WD's. For these 
WDs, the parameter hu p /kT exceeds unity long before crystal- 
lization sets in (Chabrier, Ashcroft & DeWitt 1992; Chabrier 
1993) and might well modify the diffusion coefficient by order 
unity. 
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The 22 Ne abundance in a Type la progenitor is also critical 
to the overproduction of the neutron rich isotopes 54 Fe and 58 Ni 
(Thielemann, Nomoto & Yokoi 1986). The neutron excess in 
the inner parts of an exploding WD is fixed by electron cap- 
tures. Modifying the production of these elements at that lo- 
cation would require a substantially enhanced 22 Ne abundance 
there (e.g. a pure 22 Ne core would set the initial Y e « 0.46). 
Outside of (0.3 - O.4)M , the production of 54 Fe and 58 Ni is 
fixed by the local 22 Ne abundance (Iwamoto et al. 1999). 
Hence, the production of these elements can be reduced if the 
22 Ne sinks away from the outer regions (Bravo et al. 1992). 2 
Our calculation suggests there is time for 22 Ne settling in the 
outer liquid layers of a massive WD if accreting at a low rate. 
The final answer to this puzzle awaits a time-dependent evolu- 
tion on a cosmological timescale of the diffusion equation with 
a realistic starting abundance profile and accretion history. 

There is also a 22 Ne settling calculation relevant to classi- 
cal novae. Many lines of evidence (ejected masses and overly 
abundant heavy elements) point to the ejected material in clas- 
sical novae containing matter dredged up from the underly- 
ing C/O WD (see Gehrz et al. 1998 and Starrfield 1999 for 
overviews). Livio & Truran (1994) noted that the Neon abun- 
dances measured in many of these ejecta are consistent with 
the X22 ~ 0.02 expected in a C/O WD. Hence, there is no need 
to invoke the presence of an O/Ne/Mg WD for those events. 
However, this will remain true only if the 22 Ne has not gravita- 
tionally settled out of the outer C/O layers of the WD that are 
about to be dredged up. 

We can assess this by finding how far the 22 Ne falls in the 
surface C/O layers, where g is constant. This is a simple case, 
where we integrate equation (Q) for an isothermal layer. Us- 
ing pressure as a coordinate, we find that the time it takes for 
22 Ne to fall to a pressure P starting from the surface is 



Af w 0.8 Gyr rp /3 ( - 



10 s cm s" 2 



10 20 



erg cm 
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where Tp is the value of T at pressure P and exceeds one at the 
time of ignition. 

Since the accreted material appears to dredge up a WD layer 
of mass comparable to the accreted layer, we just need to com- 
pare the 22 Ne settling time to the time it takes to accumulate the 
fuel, t a = AitGMP / g 2 M. The ratio of these timescales gives 
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2 x 10 erg cm 
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(7) 

where we have written the pressure in units of the typical igni- 
tion pressure (Livio 1994). Thus, forM > lO^'Moyr -1 (which 
applies to most classical novae, Livio 1994), the underlying 
C/O material should have 22 Ne present at X22 ~ 0.02, allow- 
ing for neon enrichment on a C/O WD in those classical novae 
that are constantly undergoing dredge-up. 

4. LUMINOSITIES FROM FALLING 22 Ne IN COOLING WHITE 
DWARFS: SEDIMENTARS 

The rate of energy release from the falling 22 Ne is given by 
the integral of the power, FV, 



FVn Ne 4irrdr. 



(8) 



Since the 22 Ne settling time is longer than the time it takes to 
cool to the onset of crystallization in the core (about 2(0.3) 
Gyr for a 0.6(1.2)M© WD; Benvenuto & Althaus 1999), we 
presume that the abundances are nearly unchanged from the 
initial state. This integral is evaluated with X22 = 0.02 and 
the core temperature, T c , being constant. For a pure oxygen 
WD at T c = 10 8 K, this gives L s =l.lx 10 29 ,3.9 x 10 29 , 1.2 x 
10 30 &4.6 x lO^ergs" 1 forM = 0.6,0.8, 1.0, & 1.2M Q . 

For massive WD's, L g becomes comparable to the cooling lu- 
minosity, L, for the same T c . This raises the important question 
of how much the settling luminosity can affect the WD cooling 
curves and thus the inferred age of the galactic disk from the 
faint end of the WD luminosity function (see Leggett, Ruiz and 
Bergeron 1998 for a recent overview). 

We begin with the common 0.6M Q WDs. The hatched re- 
gion in the top panel of Figure 2 shows L g for X22 = 0.02 as 
a function of T c . The upper (lower) bound of the hatched re- 
gion is for a pure carbon (oxygen) WD and denotes the range 
of possibilities for arbitrary C/O mixtures. The sharp downturn 
at T c w (3-4) x 10 6 K is due to core crystallization. We only 
allow for sedimentation in the liquid parts of the WD, so that, 
as T c decreases and the solid core increases, L g decreases. If 
we allowed for continuing sedimentation throughout the solid- 
ifying WD, the hatched region would just continue to the left 
with the same slope. The dashed line is the L- T c relation for a 
cooling 0.6M Q WD (see caption). For O.6M , L g is always at 
least a factor of 10 lower than the exiting luminosity, L, at that 
T c . 

However, this competition gets closer as the WD mass in- 
creases. The bottom panel of Figure 2 shows the same curves 
for 1.2M Q . In this case, L g is a factor of 5 or so less than L, 
which means that the WD cooling trajectory would be affected 
considerably. It is a coincidence that L g drops off before be- 
coming dominant and then remains nearly parallel to the L-T c 
relation. This would not occur in a pure helium WD, as crys- 
tallization never sets in. However, in that case, there would be 
no 22 Ne present. It would just be the settling of the other, less 
abundant, neutron rich metals, such as 56 Fe. Though it appears 
unlikely that the settling of iron will modify the cooling of a 
low-mass helium (M < 0.4M Q ) WD, substantial settling will 
occur in 10 Gyr. 

The point we wish to make here is that the uncertain micro- 
physics of the diffusion coefficient clearly allows for an alter- 
native view of cooling for massive WD's. Namely, it is possible 
that the settling luminosity can heat a white dwarf well enough 
to keep it in the liquid state for a very long time. For exam- 
ple, increasing D by a factor of 6 for the 1.2M Q WD would 
shift the L g curves up enough so that the WD could be powered 
by sedimentation alone: hence, our dubbing of these solutions 
as "Sedimentars". Or, if we allow for continued sedimenta- 
tion once the core is supposedly solid (say, it is supercooled), 
the lines would intersect at L m 2 x 10 30 erg s" 1 for a 1.2M Q 
WD, which would then take a time in excess of 10 Gyr to ra- 
diate all of the sedimentation energy. This sedimentar would 
live for many Gyr's at an effective temperature of 11,350 K. 
For a pure oxygen 1.36M Q WD this crossover would occur at 
L « 10 31 erg s" 1 or T e s w 21,000K, still too cold to help explain 
the possible excess of hot, massive WD's reported by Vennes 
(1999). 



2 H6flich et al. (1998, 2000) studied the differences in Type la spectral evolution that are possible if the 54 Fe was absent in the outer layers. These can be large in the 
UV, where the iron is a major opacity source. Lentz et al. (2000) varied the 54 Fe abundance in the partially burned outer layer of model W7 to see if the observations 
can accomodate a varying value. Though they also found spectral differences in the UV, nothing large enough to allow them to constrain the 22 Ne abundance. 
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5. CONCLUSIONS 

Our initial calculations suggest a number of important prob- 
lems to address in the future. The coincidence of the settling 
time with the Hubble time forces us to reconsider the physics 
of the inter-diffusion coefficients, especially in the quantum liq- 
uid realm. This is critically important for the cooling of massive 
WD's where there is some chance for the WD to shine for Gyrs 
on the power of 22 Ne sedimentation. These gravitationally pow- 
ered "sedimentars" would be a new astrophysical object that 
behave differently than the conventional cooling WD. 

The time-dependent evolution of the 22 Ne abundance in dif- 
ferent astrophysical settings is the subject of our future work. 
In that paper, we evolve both cooling and accreting WD's with 
active 22 Ne settling. This will allow for an accurate calculation 
of the evolution of the 22 Ne abundance from starting conditions 
that reflect the prior stellar evolution. For cooling WD's, we can 
follow the 22 Ne throughout the star as it crystallizes, and assess 
the possibility of altering the nature of the phase separation at 
freezing via enhanced 22 Ne abundances. 



The downward diffusion speed is very slow and it is natural 
to question whether such slow speeds are actually relevant. We 
believe so, as there is no expectation for active convection to 
undo it. Comparably weak diffusive settling also appears in the 
Sun, where inclusion of the slow gravitational settling of helium 
dramatically improves the agreement of the solar interior model 
with that measured by seismological inversions (Christensen- 
Dalsgaard et al. 1996; Basu, Pinsonneault & Bahcall 2000). 

For conversations about this problem, we thank N. Ashcroft, 
D. Chernoff, S. Kulkarni, J. Liebert, E. Salpeter, and D. Steven- 
son. For discussions about current white-dwarf cooling mod- 
els, we thank B. Hansen and D. Townsley. O. Straniero pro- 
vided information on the 22 Ne distribution in the WD's and 
the anonymous referee's comments greatly improved our dis- 
cussion. This research was partially supported by NASA via 
grant NAG 5-8658 and by the National Science Foundation un- 
der Grants PHY94-07194, PHY99-07949 and AY97-31632. L. 
B. is a Cottrell Scholar of the Research Corporation. 
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FIG. 1. — The gravitational settling time for a 22 Ne nucleus in a pure C (dashed) or pure O (solid) WD at a temperature of 10 8 K appropriate for a Type la 
progenitor. The lines show how long it takes a 22 Ne nucleus to fall to a given mass coordinate if it begins at the stellar surface. The set that begins at 1.2 (0.6) Mq 
is for a WD of that mass. 
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FIG. 2. — The luminosity from Ne settling in a 0.6Mq (top panel) and 1.2Mq (bottom panel) WD with X22 = 0.02. The hatched region shows the luminosity 
from 22 Ne settling as a function of WD core temperature. The upper (lower) bound is for a pure carbon (oxygen) WD and should bracket the range of possible C/O 
mixtures. The dashed lines are the T c —L relations for cooling DA WD's from Benvenuto & Althaus (1999). Both have helium layers of = IQ~ 2 Mq, while the 
O.6(1.2)M WD has a hydrogen layer mass of 10 _4 (10" 6 )M Q . The O.6M model is close to that of Hansen (1999) and Chabrier et al. (2000). 



